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ABSTRACT 


The  specification  of  the  surface  heat  flux  is  essential 
for  synoptic  and  seasonal  prediction  of  the  upper  ccsan 
thermal  structure.  Estimates  of  tlia  surface  heat  flux  have 
been  prepared  for  the  central  North  Pacific  during  January 
1976  through  April  1979  using  archived  fields  from  the  Fiast 
Numerical  Oceanography  Center  (FNOC)  hemispheric  atmospheric 
prediction  model.  Monthly  accumulations  of  the  surface  heat 
flux  are  compared  with  the  change  in  heat  content  above  2D0 
m  derived  from  temperature  analyses  of  the  North  Pacific  Ex- 
periment TRANSPAC  ship-of-cpportunit y  program.  Systematic 
differences  are  found  between  the  accumulated  heat  flux 
fields  and  the  oceanic  heat  content  change.  Some  of  the 
differences  are  due  to  excessively  large  changes  in  ocean 
heat  contents  above  a  fixed  level.  However,  our  earlier 
studies  have  suggested  a  bias  of  excessive  upward  surface 
heat  flux,  especially  along  the  southern  boundary  of  the  do- 
main. Assuming  local  heat  balance  o/er  a  36-month  period,  a 
correction  field  to  the  FNOC  surface  heat  flux  estimates  is 
derived.  Separate  correction  fields  for  the  heating  and 
cooling  seasons  demonstrate  a  seasonal  variation  in  the  ac- 
cumulated heat  flux  versus  heat  content  change  values. 
Thus,  six  bi-monthly  correction  fields  to  be  added  to  the 
FNOC  heat  fluxes  are  prepared  to  enable  these  heat  fluxes  to 
be    used    for   ocean    prediction. 


1.        INTRODUCTION 

Prediction  of  the  upper  ocean  the rmaL  structure  re- 
quires a  specification  of  the  surface  heat  fluxes.  The  time 
scale  on  which  the  sensible  and  latent  fluxes,  plus  the  in- 
coming and  outgoing  radiative  fluxes,  must  be  specified  is 
dependent  on  the  oceanic  phenomena  and  the  type  of  modal. 
Heat  flux  values  resolved  en  a  time  scale  of  3  hours  or  isss 
must  be  provided  if  the  iiurnal  variation  in  the  atmospheric 
forcing  is  an  important  consideration  in  the  evolution  of 
the  upper  ocaan  thermal  structure  (Garwood,  1977).  Consid- 
erable evidence  has  been  accumulated  (e.g.,  Elsberry  and 
Camp,  1973;  Camp  and  Elsberry,  197  8;  Elsberry  and  Raney, 
1978;  Elsberry  and  Garwood,  1978;  and  others)  that  the  upper 
ocean  thermal  structure  responds  significantly  to  atmospher- 
ic  synoptic-scale  forcing.  On  seasonal  time  scales,  the 
surface  heat  flux  is  accumulated  in  the  upper  ocean  layers 
during  spring  and  summer,  and  subsequently  removed  during 
autumn    and  winter. 

Ocean  mixed  layer  models  provide  a  means  of  demonstrat- 
ing the  two  ways  in  which  the  surface  heat  fluxes  affect  the 
prediction  of  ocean  thermal  structure  profiles  (Niiier  and 
Kraus,  1977).        Consider      a    well-mixed      layer    of      variable 

depth,  h,  with  temperature,  T.  The  heat  content  (H)  per 
unit   area    in   a    layer    cf   depth     (h)     is 

H   =  J    yCpT(z)dz  (1) 

where  v  is  the  density,  Cp  is  the  specific  heat,  and  z  is 
the  depth  below  the  surface.  The  change  in  heat  content  of 
the    mixed    layer   is 

3     (FCpTh)     m    ,pCfhdT   +    £CpTdh 

dt  dt  dt      ' 

because  the  changes  in  p  and  Cp  are  relatively  small.  Solv- 
ing for  the  temperature  change  from  the  first  term  on  the 
right   side   of     (2)     gives 

dT         1    d(Th)  ?    dh  1  dH 

dt         h    dt  h    dt  dC?  Ei    dt 

In  mixed  layer  models  (e.g.,  Garwood,  1977),  the  last  term 
is    written    as 
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where  n'T'H)       is      the    vertical  turbulent    heat       flux   at   the 


,,    the    Heaviside    function,    is 


base  of    the  mixed   lay  er   and  _^ 

defined    as 

dh 
_yV   =  0   for  —  <  0      ,  (4) 

dh 
=    1    for  —   >    0      . 

dt 

The  vertical  current  speed,  W  (-h)  ,  at  the  base  of  the  mixed 
layer  will   be    neglected.  This  relationship   indicates   that 

entrainment  mixing,  and  thus  the  downward  heat  flux  at  base 
of  the  layer,  is  only  associated  with  deepening  layers.  The 
fraction  of  turbulent  kinetic  energy  that  is  available  for 
entrainment  aixing  in  the  Garwood  model  is  dependent  on  both 
the  surface  friction  velocity  and  the  surface  buoyancy  flux, 
which  is  determined  partly  by  the  sjc  face  heat  flux.  Conse- 
guently,  the  surface  heat  flux  contributes  to  both  terms  in 
(3).  The  effect  in  the  first  term  is  reflected  directly  in 
the  change  in  heat  content  due  to  the  surface  flux.  The 
thermal  structure  is  strongly  dependent  on  the  vertical  re- 
distribution of  the  heat  via  entrainment  mixing,  which  is 
partly    due   to    the    upward   surface  heat    flux. 

We  conclude  that  specification  of  the  surface  heat  flux 
is  an  essential  factor  for  prediction  of  anomalous  ocean 
thermal  structure.  In  the  next  section,  we  wiLi  briefly  re- 
view the  methods  available  for  estimating  the  surface  hea-1- 
flux  over  the  ocean.  In  the  following  section,  we  examine 
the  role  of  the  surface  heat  flux  in  the  oceanic  heat  budg- 
et. We  then  derive  a  correction  field  to  renove  a  bias  in 
the  surface  heat  flux  that  would  be  detrimental  to  ocean- 
thermal    structure  prediction. 

2.       SPECIFICATION    OF    THE    SH8 FACE    HEAT    FLUX 
The    bulk    aerodynamic    method      for    calculating    the    sensi- 
ble and    latent      heat    fluxes  involves    the      surface    wind    field 
and  the     iifferences    in    temperature    and    specific    h;imidity    be- 
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tween  the  sea  and  the  air.  Calculations  of  tha  incoming  and 
outgoing  radiative  fluxes  reguire  knowledge  of  the  cloud 
cover  as  well  as  the  solar  altitude.  An  exampla  of  a  system 
for  estimating  the  surface  heat  flax  from  ship  observations 
is    given    by   Clark,    et    al.  (1974),    Bunker    (1976)       and   Clark 

(1981).  A  common  approach  is  to  calculate  tha  surface  flux 
from  each  ship  report  and  then  average  over  some  space  and 
time  interval.  For   example,      all   ship      observations   taken 

during  a  month  within  5  latitude  and  longitude  may  be  used 
to  represent  the  average  heat  flux  in  that  domain.  Husby 
(1980)  has  compared  the  anomalous  h  ear:  flux  over  six-month 
periods  at  Ocean  Weather  Station  V  (3U*N,  164°3)  with  the 
estimate  based  on  merchant  ship  observations  within  a  4fl 
guadrangle.  He  concludes  that  the  merchant  ship  reports 
must  be  carefully  screened  prior  to  the  heat  flux  computa- 
tions. 

The  accuracy  of  this  method  is  clearly  dependent  on  the 
number  of  observations.  There  is  no  assurance  that  the  re- 
ports will  be  randomly  distributed  in  space  or  time.  One 
problem  appears  to  be  a  "fair  weather  bias"  because  the 
ships  tend  to  avoid  bad  weather.  3arr.ett  (1931)  indicates 
that  a  set  of  heat  fluxes  in  the  central  North  Pacific  Ocean 
estimated  by  this  method  had  a  bias  of  30-U5  H  m"  ,  which 
would  result  in  an  excessive  estimation  of  the  heat  flux 
into   the    ocean. 

An  estimate  of  the  monthly  heat  flux  is  not  sufficient 
for  short-term  ocean  prediction,  which  reguirss  information 
on  diurnal,  or  at  least  synoptic,  time  scalas.  Synoptic 
maps  of  heat  flux  are  difficult  to  analyze  because  of  the 
sparsity  of  ship  observations.  An  indirect  method  is  pro- 
posed here.  The  atmospheric  prediction  models  also  require 
a  calculation  of  the  heat  flux  at  the  ocean  surface.  The 
heating  package  for  the  Fleet  Numerical  Oceanography  Center 
(FNOC)  model  involves  bulk  aerodynamic  calculations  of  the 
latent  and  sensible  heat  fluxes,  plus  radiati7e  fluxes  that 
are  a  function  of  the  model-estimated  cloudinsss  (Kesei  and 
*  inning h off ,     1972). 
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Our  basic  hypothesis  is  that  the  FNOC  atmospheric  pre- 
diction model  heating  package  can  provide  the  heat  fluxss 
necessary      for      ocean      prediction     experiments.  Gallacher 

(1979)  has  described  the  method  used  for  extracting  hourly 
heat  flux  estimates  from  the  FNOC  archives.  Elsberry,  Gal- 
lacher and  Garwood  (1979)  have  used  these  heat  fluxes  to 
predict  the  ocean  thermal  structure  changes  during  the  au- 
tumn of  1976.  3udd  (1980)  also  used  these  fields  in  an  at- 
tempt to  predict  the  spring  transition  from  the  winter  to 
the  summer  regimes  in  the  central  North  Pacific  Ocean.  Budd 
found  a  systematic  bias  with  a  too  large  upward  heat  flux 
near  30*N.  Recently,  Steiner  (198  1)  also  found  a  systematic 
bias  in  these  hear  flux  estimates  in  the  region  between  Ha- 
waii and    San   Francisco. 

The  purpose  of  this  paper  is  to  describe  a  correction 
field  to  the  FNOC  surface  heat  flux  fields  to  permit  their 
use   in    ocean    prediction   experiments.       As   in   the    case    of   Budd 

(1980)  and  Steiner  (1981)  ,  the  approach  is  to  derive  a  cor- 
rection field  that  assures  a  long  term  (annual  or  longer) 
heat  balance  in  the  upper  ocean.  Analyses  by  White  and  Bern- 
stein (1979)  of  the  ocean  thermal  structure  observations  in 
the  TRANSPAC  ship-of- opportunity  program  are  used  to  calcu- 
late the  time  changes  in  oceanic  heat  content.  These  are 
compared  with  the  accumulated  surface  heat  fluxes  derived 
from  the  FNOC  files.  We  then  derive  the  correction  to  the 
monthly  surface  heat  flux  that  is  necessary  to  assure  local 
heat  balance  over  the  1976-1978  period.  In  the  next  section 
we  review  briefly  some  tests  of  the  local  heat  balance  as- 
sumption. It  is  only  on  the  very  long  time  interval  that  we 
assume  local  heat  balance.  We  do  not  reguire  or  insure  lo- 
cal hear  balance  en  monthly  or  shorter  time  intervals  over 
which   we    are    doing   prediction    experiments. 


3-       OCEAN    HEAT    BUDGET    STUDIES 

The    general    purpose    in  this  section   is   to    determine   the 
relative    importance      of   the  surface      heat    flux    in      the    upper 


ocean   he  at    budcet. 


Our 


'imary    interest    is 


coen-ccean 
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regimes  that  would  be  similar  to  tha  conditions  in  the  North 
Pacific  Experiment  (NORPAX)  Anomaly  Dynamics  Study  (ADS)  re- 
gion. We  do  no+  consider  regions  of  boundary  currents  or 
near-equatorial  areas  in  which  horizontal  and  vertical  ad- 
vection  are  likely  to  be  significant.  The  space  and  time 
scales    of    interest    are   greater   than    1  000    km   and    one    month. 

Many  comparisons  have  been  made  betwean  the  seasonal 
changes  of  heat  content  of  the  upper  ocaan  and  the  accumu- 
lated heat  flux  at  the  surface.  Bryan  and  Schroeder  (1950) 
compared  the  heat  content  calculated  from  North  Atlantic  BT 
data  with  the  surface  heating  estimated  by  Budyko  (1955). 
They  found  that  the  surface  heating  on  a  seasonal  basis  was 
about  20%  less  thar  the  change  in  halt  content  in  the  region 
between  20°N  and  about  50°N.  By  contrast.  Bather.  (1971) 
found  that  the  surface  heat  exchange  estimates  of  Wyrtki 
(1966)  could    account      for  only      29$    of      tha    local      monthly 

change  in  heat  content  in  the  North  Pacific  Dcean.  Gill  and 
Niiler    (1973)  propose    that      the    heat     input    averaged      over 

large  areas  and  times  is  mainly  stored  locally,  and  horizon- 
tal advection  by  the  mean  flow  is  not  particularly  impor- 
tant. They  also  cite  comparisons  by  Tabata  (1965)  and  Rob- 
inson (1966)  using  Ocean  Weather  Ship  data  which  suggest 
that  most  of  the  heat  input  changes  ire  stored  locally. 
Gill  and  Niiler  further  suggest  that  the  inaccuracy  of  the 
heat  flux  estimates  over  large  areas  away  from  Ocean  Weather 
Ships  is  the  likely  cause  of  some  of  the  departure  from  lo- 
cal  heat    balance. 

A  careful  study  of  the  upper  ocean  heat  budget  near  OWS 
?  based  on  two  weeks  of  high-quality  observations  during  the 
Mixed  Layer  Experiment  (MILE)  has  bean  re  portal  by  Davis,  et 
ai.  (198  1).  They  found  that  a  one- dimensional  upper  layer 
heat  budget  may  be  closed  acceptably  if  the  temperature  is 
well-sampled,  and  if  the  vertical  velocity  in  the  seasonal 
thermocline   is    also   taken    into    account. 

A    recent    study    by    3arnett     {1 981  •        used    AX3T's    along    158° 
and    170*W    between    30*    and      50° N   at    approximately    monthly    in- 
tervals   to   estimate   the   heat    budget.       Barnett    concludes   that 
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approximately  90-95%  of  the  variance  in  the  seasonal  change 
of  heat  storage  can  be  accounted  for  by  the  air/sea  heat  ex- 
change. Horizontal  and  vertical  advection  were  of  limited 
and  no  use,  respectively,  in  reproducing  the  seasonal  cycle 
variance.  Barnett  also  concludes  that  prior  heat-budget 
studies  that  used  an  arbitrary  lower  depth  in  the  estimation 
of  heat  content  may  have  inadvertently  included  variance  on 
different  space  and  time  scales.  Finally,  Barnett  suggests 
that  the  inability  to  explain  adequately  the  anomalous  heat 
content  changes  in  the  central  Pacific  is  likely  due  to  im- 
precise   knowledge    cf   the   source   terms    in   the    heat    budget. 

Oar  earlier  studies  (Elsberry  et  al.  ,  1979;  Budd,  1989; 
Steiner,  1981)  have  used  an  oceanic  heat  budget  to  determine 
the  feasibility  of  using  the  FNOC  heat  flux  estimates  for 
ocean  prediction.  In  each  case,  the  imbalances  appear  to  be 
systematic  in  space  and  time.  Large,  short-term  imbalances 
are  likely  due  to  ocean  observational  sampling  errors  and  to 
using  an  arbitrary  lower  depth  in  the  estimate  of  heat  con- 
tent (Barnett,  1981).  However,  the  long-term,  systematic 
error  is  likely  due  to  a  bias  in  the  FNOC  heat  flux  esti- 
mates. In  the  following  sections,  *e  derive  an  appropriate 
correction  field  to  be  added  to  the  FNOC  heat  flux  estimates 
to    remove   the    long    term   bias. 

4.  CALCULATION  CF  THE  0C2AN  HBAT  BUDGET 
Monthly  mean  temperatures  during  1976-1979  have  been 
objectively  analyzed  by  White  and  Bernstein  (1979)  based  on 
TRANSPAC  ship-of-opportuni ty  XBT's  only.  Analyzed  tempera- 
tures are  available  on  a  2°  latitude  by  5*  longitude  grid  at 
0,  20,  U0,  60,  90,  120,  15  0,  200,  250,  300  ani  U00  m.  The 
general  domain  cf  the  heat  budget  calculations  is  from  30* 
to  50*N  and  from  170°S  to  130°W.  However,  the  number  of  ob- 
servations in  the  northwest  and  southwest  corners  of  the 
grid  are  insufficient  to  provide  reliable  estimates  during 
some  months.  Consequently,  these  areas  are  eliminated  from 
the  following  analysis,  and  will  appear  as  cross-hatched 
areas   on    all    figures. 
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To  calculate  the  heat  content  as  in  (1),  the  objective- 
ly analyzed  temperatures  are  first  interpolated  to  5  m 
depths.  The  heat  content  is  actually  calculated  relative  to 
the  200  D  temperature,  with  the  intent  of  renoving  some  of 
the  change  in  heat  content  due  to  vertical  displacements  of 
the  thermccline.  Barnett  (1981)  has  suggested  that  selec- 
tion of  an  arbitrary  depth  such  as  2D0  m  may  introduce  addi- 
tional variaice  into  the  heat  content  calculation  due  to  in- 
ternal waves,  tides  and  so  forth.  Inspection  of 
month-to-month  heat  content  values  does  reveal  time  changes 
which  are  much  larger  than  can  be  reasonably  accounted  for 
by  surface  heat  fluxes.  Large  positive  time  changes  tend  to 
be  followed  by  large  negative  changes.  Thus,  time  intervals 
of  either  two  or  three  months  have  been  teste!  in  the  esti- 
mate of    the  oceanic   heat   content  change, 

£H    =    H(t    ♦    ^t)      -   H(t)     .  (5) 

We  expect  ,/^H  to  be  positive  during  the  period  of  net  down- 
ward surface  heat  flux  (roughly  between  April  and  September) 
and  negative  during  the  period  of  net  upward  flux.  The 
monthly  TRANSPAC  analyses  are  assumed  to  apply  on  the  15th 
of  each  month.  If  ^t  is  two  months  in  (5),  then  /\H  will 
also  apply  at  the  15th  of  the  month.  For  example,  the  dif- 
ference between  the  March  1976  and  the  January  1976  heat 
content  fields  would  be  applied  at  the  middle  of  February. 
Notice  that  this  is  the  first  possible  value  that  can  be 
calculated  because  the  analyses  begin  in  January  1976.  How- 
ever, if  £^t  is  three  months  in  (5),  than  ^H  will  apply  at 
the   mid-pcint    of   this    interval.  Thus,      the    first    possible 

difference  would  be  between  April  1976  and  January  1976,  and 
it  would  apply  on  1  March  1976.  The  advantage  of  the  longer 
time  interval  is  that  shorter  period  fluctuations  in  heat 
content  tend  to  be  averaged  out  and  one  obtains  a  more  con- 
sistent measure  of  the  seasonal  variation.  The  disadvantage 
of  the  longer  time  interval  is  the  greater  inaccuracy  in  the 
finite  difference  approximation  to  a  derivative.  Finally, 
it  is  necessary  to  extend  the  ^H  calculation  into  1979  to 
obtain    three    complete    annual    cycles.       In   the    case    of    £*-    =    2 
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(3)  months,  the  analyses  through  February  (March)  1979  are 
used. 

Following    Emery    (1976),    the  basic  hsat   balance    equation 
may  be    written    as 

Be      wd  (h  -  dcpdt  ) 

--  =    -£  *  p  -    V  •  *CH    ♦   Q    ,  (6) 

dt         D 

where  H  is  defined  as  in  (1)  with  -h ,  the  mixed  layer 
depth,  replaced  by  D  which  is  a  (constant)  depth  below  the 
main  thermocline.  In  (6 )  ,  the  first  term  on  the  right  rep- 
resents the  divergence  of  heat  due  to  horizontal  divergence, 
wpH/D,  and  vertical  advection  W0pCpT  .  The  second  term  is 
the  horizontal  advection  of  heat  and  the  last  term  is  the 
surface    heat   flux.  The    suggestion    by   Barnett       (1981)       and 

Davis  et .  al.  (1981)  that  D  be  the  depth  of  an  isotherm  be- 
low the  main  thermocline  rather  than  a  constant  depth  re- 
sults in  the  divergence  term  in  (6)  being  z?ro,  and  thus 
eliminates  a  source  of  uncertainty  from  the  heat  budget. 
Due  to  the  large  vertical  sampling  interval  in  the  TRANSPAC 
analyses,  we  will  use  D  =  200m  throughout  this  study.  Thus, 
we    assume  the    following   budget    equation 

/"■At 

&H   =   J    Qdt    +    Residual  =    QT   ♦   Residual    .  (7) 

The  integral  sign  indicates  that  ths  air-sea  flux  is  summed 
over  the  same  time  interval  that  &$  is  evaluated.  QT  is 
defined  as  the  surface  heating.  In  the  ^t  =  2  month  exam- 
pie  above,  the  sum  would  be  from  15  January  through  15 
March.  The  residual  term  in  (7)  includes  the  remaining 
terms  (physical  effects)  in  (6)  that  can  not  be  evaluated, 
plus  the  errors  in  estimating  the  haat  content  changes  and 
the  surface  fluxes.  If  these  observational  errors  are  ran- 
dom over  a  sufficiently  long  interval,  their  contribution 
should    be      averaged    cut.  In    partioular,         w=    sum      over    36 

months  (or  a  small  subset  if  a  seasonally  varying  correction 
factor  is  derived)  and  derive  an  averaged  correction  field 
to  be  added  to  the  surface  heat  flux.  However,  ^he  physical 
terms  may  not  sum  to  zero.  Persistant  vertical  and  horizon- 
tal advection  or  diffusive  effects  that  accumulate  during 
the    period     are    thus   also      incorporated  into      khe    correction 


field.  The  impact  of  assuming  a  local  (one-dimensional) 
heat  balance  via  this  long-term  correction  fieLd  will  be  as- 
sessed  in   the    prediction   experiments. 

The    total    heat    flux   Q    in    (7)     may    be    expressed   as 

C  =   Q4    -    (Qb   ♦   Qw    +    Q1)    .  (8) 

where  the  subscripts  s,  b,  1  and  h  refer  to  solar,  back,  la- 
tent and  sensible  heat  fluxes  through  the  sea  surface.  A 
description  of  the  calculation  of  each  of  these  heat  fluxes 
in  the  FNOC  atmospheric  prediction  model  is  given  in  Kesel 
and  Winninghoff  (1972).  The  values  cf  these  terms  after  a 
one-hour  integration  of  the  model  are  taken  to  apply  at  the 
synoptic  map  time  (00  or  12  GMT).  After  spatial  interpola- 
tion to  the  grid  points  of  the  ocean  analysis,  these  instan- 
taneous values  are  interpolated  to  hourly  values  using  tech- 
nigues  described  by  Gallacher  (1979).  The  representation  of 
the  solar  flux  is  particularly  involved.  It  is  this  time 
series  of  hourly  values  d f  Q4  and  Qu  (=  Q,  f  2K  +  Q, )  that 
are  accumulated  in  (7).  Simply  summing  the  instantaneous 
12-h  FNOC  values  would  not  produce  the  same  sum  because  we 
have   reconstructed    the    hourly    variation   in    Q. 

An  example  of  a  field  of  ^H  near  the  beginning  of  the 
ocean  heating  season  is  given  in  Fig.  1  .  Dver  most  (ail) 
of  the  domain,  the  heat  content  change  is  positive,  as  ex- 
pected.         A  value      of   2    x    10      cai    cm  over    this   twc-month 

period  corresponds  to  an  average  increase  in  temperature  of 
2°  C  over  a  depth  of  100  m,  or  &°Z  over  a  depth  of  50  m. 
Less  confidence  should  be  given  to  the  values  near  the 
northern  and  southern  boundaries  of  the  domain.  The  primary 
ship-cf-opport unity  tracks  lie  on  the  great  r-ircie  route  be- 
tween Japan  and  the  USA  west  coast  ports  (White  and  Bern- 
stein, 1979).  During  the  winter,  the  tracks  tend  to  be  dis- 
placed southward,  which  makes  the  northern  region  sampling 
rather  poor.  Somewhat  the  opposite  effect  occurs  during  the 
summer  as  the  ship  tracks  are  extended  northward.  The  3b- 
iective  analysis  procedure  (White  and  Bernstein,  1979)  used 
through  1978  would  generate  values  at  all  gridpoints  regard- 
less of    whether    any    observations    fell    in    the    vicinity   of    -.he 
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point.  Because  the  deficient  points  tend  to  lie  toward  the 
north  and  the  south,  the  fictitious  values  tend  to  resemble 
northward  or  southward  extrapolations  from  the  central  re- 
gion. 

The  corresponding  integrated  surface  heat  flux  for  this 
period  derived  from  the  FNOC  calculations  is  shown  in  Fig. 
1c.  In  contrast  to  the  &E  fields,  no  extrapolation  of  the 
fields  toward  the  domain  boundaries  is  involved  in  the  Qr 
field.  In  the  band  from  40-50°  N,  the  surface  exchange  is 
approximately  egual  to  the  observed  ocean  heat  content 
change,  which  suggests  an  approximate  local  heat  balance. 
However,  the  net  FNOC  heat  flux  continues  to  be  upward  be- 
tween 30-35  N  during  this  period.  One  can  not  tell  from 
this  diagram  whether  the  daily  values  of  Q5  in  (8)  are  too 
small,  or  if  the  upward  surface  heat  flax  (Qw)  is  too  large. 
The  fact  that  the  maximum  values  are  found  near  the  longi- 
tude for  which  the  boundary  of  the  FNOC  grid  is  tangent  to 
the  eguatcr  (therefore  closest  to  the  study  region)  suggests 
that  the  problem  may  lie  in  the  boundary  conditions  that  are 
applied. 

The  difference  between  £\  H  and  the  cumulative  heat  flux 
(Q_)  is  shown  in  Fig.  1  .  A  light  filter,  with  weights 
2-4-2  in  latitude  and  1-4-1  in  longitude,  has  been  passed 
over  this  field  to  reduce  small  scale  noise.  Mt  hough  there 
are  considerable  areas  with  the  expected  zero  values,  there 
are  other  areas  with  positive  and  negative  valies.  In  par- 
ticular, the  residual  term  in  (8)  is  very  iirge  along  the 
southern  boundary  because  of  the  lack  of  a  net  downward  sur- 
face heat    flux     (Fig.       1    )  . 

The  complete  set  of  monthly  ^3  ,  2_  and  £*H-Q_  values 
is  contained  in  the  Appendix.  Since  our  desire  is  to  obtain 
a  correction  field,  it  is  not  appropriate  to  discuss  here 
each  of      these    sets.  One  general      feature    is      that    the      QT 

fields  have  less  east-west  variation  ana  nore  north-south 
variation  than  do  the  ^H  fields.  There  are  several  possi- 
ble causes  for  this  feature:  (a)  the  FNOC  heat  fluxes  may 
not   contain   an    adeouate   representation    of    the    =ast-west    var- 
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iations  in  the  anomalous  surface  flaxes;  (b)  vertical  advec- 
tive  processas  associated  with  ocean  eddies  are  contained  in 
the  heat  content  values  relative  to  a  fixed  depth;  and  (c) 
errors  due  to  inadequate  observations  or  analysis  of  the 
ocean  heat  content.  The  spring  transition  regime  example  in 
Fig.  1  illustrates  one  of  the  worst  imbalances  in  the  time 
series.  An  excessive  upward  heat  flax  also  occurs  along  the 
southern  boundary  during  autumn.  However,  the  ^  H  field  is 
negative  (upper  ocean  heat  content  diminishing  in  time)  dar- 
ing this  period.  Consequently,  ths  magnitude  and  sign  of 
the  imbalance  may  have  a  seasonal  dependence  (discussed  fur- 
ther  in    the  next   section). 

Similar  sets  of  & H,  Q-  and  A H-QT  were  also  prepared 
with  ^t  =  3  months.  These  fields  (not  shown}  are  not  di- 
rectly comparable  with  the  two-month  set  because  of  the 
half-month  displacement  of  the  central  point.  However,  the 
time  evolution  of  the  major  features  using  ^\  t  =  3  months 
can  be  easily  associated  with  the  evolution  shown  by  the  £\t 
=  2  month  series.  Therefore,  only  the  correction  fields 
based   on    the    2-month    differences   will    be   discussed    below. 

5.       INTERPRETATION    OF   THE    HEAr     BUDGET    IMBALANCE 

Perhaps  the  most  important  factor  to  emphasize  regard- 
ing these  calculations  of  the  local  heat  budget  (8)  is  that 
each  of  the  two  terms  is  subject  to  large  errors.  As  noted 
above,  calculating  the  heat  content  changes  relative  to  a 
fixed  depth  may  produce  oscillations  equivalent  to  a  1-2  C 
change  averaged  over  the  entire  depth.  A  faulty  XBT  could 
also  be  the  cause  of  a  temperature  bias  with  depth.  How- 
ever, we  would  expect  that  the  objective  analysis  technique 
would  tend  to  eliminate  such  a  bias  if  other  correct  pro- 
files are  in  that  region.  3ased  on  our  earlier  studies 
(Elsberry,  et  al.,  1979  ;  3udd,  1930;  Steiner,  1981),  we  ex- 
pect that  the  monthly  surface  heating  dees  not  have  a  suffi- 
ciently large  seasonal  amplitude.  There  is  also  a  persis- 
tent bias  toward  excessive  heat  loss  to  the  atnosphere  along 
the    southern    boundary   of   the    domain. 
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Because  the  heat  budget  imbalances  may  arise  from  ei- 
ther term,  Table  1  was  prepared  to  indicate  the  various  pos- 
sibilities that  may  arise.  It  is  useful  to  separate  the  AH 
into  separate  periods  when  the  ocean  heat  content  is  in- 
creasing cr  decreasing.  Apart  from  small  calculation  er- 
rors, one  would  hope  that  the  corresponding  QT  would  be  po- 
sitive (surface  heating)  or  negative,  respectively,  so  that 
Cases  A-B  or  E-F  would  apply.  One  axpects  approximately  an 
equal  distribution  between  A  and  B  or  between  E  and  F  if  the 
errors  are  random.  The  percentage  of  gridpoints  with  month- 
ly  ^H-QT  differences  exceeding  0.5  x  10  cal  cm  are  nor 
evenly  distributed  between  cases  A  an d  3  or  between  E  and  F 
when  the  uncorrected  QT  is  used  in  the  differences.  There 
is  clearly  a  bias  toward  cases  A  ani  E,  which  could  be  at- 
tributed to  excessive  upward  heat  flux  during  both  the  ocean 
warming  and  cooling  periods.  However  an  approximately  equal 
distribution  is  obtained  when  the  corrected  (see  description 
of  six  bi-monthly  correction  fields  in  the  next  section) 
surface    heat    flux   is   used    in   the  differences. 

A  value  of  0.5  x  10  cal  cm  corresponds  :o  a  heat  flux 
bias  of  7.0  cal  cm  hr  or  a  temperature  bias  of  .25*  C  in  a 
200  m  water  column.  Both  of  these  values  are  within  the  ex- 
pected range  of  instrument  error  for  the  measurements  which 
were  used  in  the  analyses.  Whereas  55X  of  all  the  differ- 
ences exceed  this  criterion  for  the  uncorrected  heat  flux, 
only   42^    exceed    the    criterion    for   the    corrected    heat    flux. 

Cases  C  and  D  in  Table  1  are  labelled  as  drastic  imba- 
lances because  the  ^H  and  Qr  are  of  opposite  tendencies.  A 
physical  explanation  for  such  an  event  might  be  the  advec- 
tion  of  a  warm  (cold)  ocean  eddy  into  the  region  that  has 
upward  (downward)  surface  heat  flux.  An  example  of  a  compu- 
tational explanation  is  found  along  the  southern  boundary  in 
Fig.  1  where  QT  <  0  when  ^ H  >  0.  The  percentages  in  Table 
1  indicate  that  drastic  imbalances  are  relatively  rare  dur- 
ing ocean  cooling  periods  (Case  D  -  1  %)  when  the  uncorrected 
QT  is  used  in  the  differences.  However,  this  is  not  the 
case  during  the    ocean    warming    periois,    when    a.    large    fraction 
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of  the  points  have  significant  upward  heat  fluxes  rather 
than  downward.  Since  the  heat  budget  for  the  central  Pacif- 
ic is  in  an  approximate  local  balance  for  seasonal  time 
scales  (Gill  and  Niiler,  1973),  the  percentages  for  Cases  C 
and  D  should  be  less  than  those  for  Cases  A,  3,  E  and  F. 
This  distribution  is  obtained  when  the  corrected  QT  values 
are  used,  but  no-  with  the  uncorrected  0_T.  This  is  a  fur- 
ther demonstration  of  the  need  to  adjust  the  FNOC  surface 
heat   fluxes. 

The  alternative  of  a  systematic  bias  due  to  the  AH 
calculation  does  not  appear  likely.  Although  there  are  more 
points  with  £H  positive  (59%)  than  negative  (41%),  this  as- 
ymmetry is  probably  within  the  limits  of  roughly  offsetting 
periods  of  ocean  warming  versus  cooling.  The  imbalances  be- 
tween seasons  in  Table  1  are  not  consistent  with  horizontal 
advection  being  a  primary  cause.  One  would  expect  stronger 
Skman  advection  effects  during  the  winter  (cooling  season.), 
whereas    the   larger    imbalances    are    found   during    the   summer. 

TABLE    1 
A  H  versus    QT 


(Downward   surrace 
heat    flux) 

Qt    <   0 
(Upward   surface 
heat    flux) 

AH    >    0 

Ocean    warming 

Case    A           AH-QT    >0 
Inadequate    downward    Q 
Excessive  ocean   warming 
22.4%            (26.8%) 

Oa=e    C         AH-2t     >    ° 

Drastic   imbalance 

3  3.4%           (7.3%) 

Case    B         AH-QT   <    0 
Excessive  downward    Q 
Inadequate    ocean    warming 
2.9%              (20.0%) 

rt      A^    <    0. 
Ocean    coding 

Case    D        AH-QX    <    0 

Drastic    imbalance 

1.056              (6.5%) 

Case    E         £±'1-QT    >    0 
Excessive    UDward    Q 
Inadequate    ocean    coclir. 
3  a.  6%           (16.7%) 

Case    ?        AH-3r    <    0 
Inadequate   upward   Q 
Excessive   ocean   cooling 
5.4?5              (22.6%) 

Possible    physical    association    for    positive   or    negative    heat    budget 
imbalances   given  that   the    heat    content   or    integrated    surface    fluxes 
are   oositive   or    negative.       Percentage    of   the   gridpoints    for    the    36 
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6.       CALCULATION  OF    CORRECTION    FIELD 

Our  objective  is  to  determine  a  correction  field  that 
may  be  applied  to  the  Q  values  in  (3)  during  the  ocean  pre- 
diction experiments.  Because  the  time  step  in  the  Garwood 
model   is    typically    1    h,      it  is   desirable  that   the    correction 
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be  in  units  of  cal  cm  hour  .  This  is  done  by  converting 
each  of  ^H-QT  fields  to  an  hourly  rate  using  the  proper 
number   of      hours   in      the   time    interval.  The    sign      is   also 

changed      because  the     correction  is      to      be   added      to   the      Qy- 
field. 

It  is  not  within  the  scope  of  this  paper  to  determine 
the  fraction  of  this  correction  which  will  be  applied  to  Qs 
versus  the  remaining  three  terms  in  (8).  It  night  be  noted 
that  Qs  occurs  only  during  the  daytime  and  this  flux  is  to 
be  distributed  exponentially  with  iepth  by  the  prediction 
model.  In  contrast,  the  remaining  three  terms  in  (8)  apply 
only  at  the  surface  and  tend  to  be  positive  (upward) 
throughout  the  day.  Thus  ,  one  may  expect  considerable  dif- 
ferences in  the  upper  ocean  predictions  as  the  fraction  of 
the   correction   that    is  allcted    to   Qs    is  increased. 

For  simplicity,  one  desires  a  single  correction  field 
as  in  Fig.  2.  The  basic  features  in  this  correction  field 
are  generally  consistent  with  the  pattern  shown,  in  Fig.  1a. 
In  particular,  the  correction  reduces  the  surface  heat  flax 
(upward  is  defined  to  be  positive)  along  the  southern  bound- 
ary. It  can  be  shown  that  this  correction  is  equivalent  to 
calculating  the  residual  in  (7)  using  the  difference  in  H 
between  January-February  1979  and  January-February  1976,  and 
the  Q_  for  the  entire  period.  We  nay  also  regard  the  single 
correction  field  as  the  adjustment  necessary  during  each 
hourly  time  step  of  an  integration  from  15  January  1976  to 
15  January  1979.  If  the  correction  is  applied  in  this  way, 
we  can  be  assured  of  conservation  of  heat  at  =ach  gridpoint 
for  the  entire  three  year  period  (assuming  no  vertical  dif- 
fusion  and   no    round-off  errors). 

Each  of  the  36  naps  of  AH~'2r  -s  ^  — -2-5n~-  from  that 
implied    by  the    single   correction    field.      We    noted    above   that 
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the  seasonal  variation  in  ^\H  or  QT  a ight  be  associated  with 
a  modulation  in  the  difference  field.  A  separate  correction 
field  is  derived  for  November  through  April  (Fig.  3a)  and 
for  May  through  October  (Fig.  3b).  Each  of  these  correc- 
tions is  summed  over  the  three  years  and  six  months,  so  that 
a  total  of  18  values  is  included  in  each  average.  There  is 
clearly  a  seasonal  variation  in  these  two  correction  fields. 
The  pattern  during  the  heating  season  (Fig.  3b)  is  similar 
to  the  single  correction  field  (Fig.  2),  but  the  values  are 
larger.  The  pattern  during  the  cooling  season  (Fig.  3a)  is 
quite  different.  Not  only  are  the  values  considerably 
smaller  than  in  Fig.  3br  there  is  mere  of  a  north-south  ori- 
entation   of  the    isclines    in   the   north-central    region. 

Based  on  the  seasonality  in  Fig.  3,  we  also  examined 
further  subdivisions  into  quarterly  or  bi-monthly  correction 
fields.  A  decision  was  then  made  to  adopt  the  six  correc- 
tion fields  in  Fig.  U  a- f .  It  is  felt  that  six  maps  will 
give  a  better  representation  of  the  seasonal  variation.  In- 
spection of  Fig.  4  b-e  reveals  a  basically  east-west  pat- 
tern, whereas  the  remaining  two  maps  exhibit  the  north-south 
orientation.  Even  though  the  seasonal  variation  does  not 
appear  to  be  a  sine  wave,  the  six  maps  appear  to  provide  a 
relatively  smooth  transition  between  the  two  basic  patterns 
in  Fia.  3.  There  are  some  non-seasonal  features  in  these 
correction  fields.  The  largest  is  found  along  the  western 
side  of  the  domain.  Larger  corrections  are  required  during 
February-March  (Fig.  4  a)  and  during  August-September  (Fig. 
4 d) .  The  transition  between  the  correction  field  for  Decem- 
ber-January (Fig.  4e)  and  for  Febri ary-Maroh  is  especially 
noteworthy  along  the  western  boundary.  The  regions  of  posi- 
tive corrections  (additional  upward  heat  flux  required)  also 
tend  to  be  somewhat  erratic.  Howavar,  these  positive  cor- 
rections are  always  small,  so  they  will  have  little  effect 
on    the    ocean   predictions. 

An  example  cf  the  application  3f  the  heat  flux  correc- 
tion field  at  a  specific  location  is  given  in  Fig.  5.  The 
seascr.al      changes      in      heat        content       calculated      from      the 
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TRANSPAC    analyses   have  relatively   constant   values    during  the 
heating    season.  By  contrast,      ths      ^H    values      during   the 

cooling  season  tend  to  reach  rather  well-defined  peaks.  The 
surface  heating  agrees  rather  well  with  the  maximum  ^  H  val- 
ues during  each  coding  season.  However,  the  uncorrected 
heat  flux  during  the  warming  season  is  clearly  deficient, 
especially  during  1978.  Applying  the  heat  flax  corrections 
for  this  location  from  Fig.  4  improves  the  agreement  between 
the  surface  heating  and  the  ^  H.  The  major  feature  to  be 
noticed  is  that  the  (^ H  -  QT)  differences  are  indeed  sys- 
tematic, so  that  a  single  correction  in  each  two-month  peri- 
od tends  to  improve  all  three  years.  One  finds  periods  in 
which  the  ^H  -  QT  remains  large  after  the  correction  has 
been  applied.  This  indicates  that  although  we  have  had  to 
correct  the  surface  heat  flux  we  have  not  forced  the  heat 
budget  to  be  one-dimensional  for  periods  shorter  than  36 
months.  Some  examples  are  the  differences  daring  April  - 
June  1976  and  during  July  -  August  1978.  Another  feature  in 
Fig.  5  is  that  the  area  between  the  ^H  and  the  corrected  QT 
curves  must  sum  to  the  difference  between  the  final  and  ini- 
tial H  values  (February  1979  and  January  1976,  respective- 
ly). This  reguirement  is  a  consequence  of  the  local  heat 
budget    assumption    over   the   36-month   period. 

7.       SUMMARY 

We  have  prepared  six  bi-monthly  correction  fields  to  be 
applied  to  the  FNOC  heat  flux  valuas  to  be  used  for  ocean 
prediction.  The  largest  corrections  are  found  generally  be- 
tween 30  and  38  N  during  April  through  November.  The  heat 
flux  bias  is  evidently  not  serious  for  the  atmospheric  pre- 
dictions because  they  are  limited  to  72  h.  However,  such  a 
bias  can  be  disastrous  (Budd,  1980>  for  ocean  prediction 
over  monthly  time  periods.  The  correction  fields  are  aver- 
aged over  three  years  (1976-1979).  Closure  of  the  local 
heat  budget  over  the  entire  period  is  insurad  by  the  use  of 
this  correction  field.  However,  this  is  not  true  over  any 
shorter    time    intervals.       I?,   particular,    fluctuations    in    heat 
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content    in  the      monthly  TRANSPAC  analyses    are      not   accounted 
for  by   the  long-term   correction   field. 

The  success  of  these  correction  fields  can  only  be 
judged  by  their  application  in  the  ocean  prediction  experi- 
ments. Inclusion  of  the  corrections  should  result  in  im- 
proved predictions.  One  of  the  tests  will  be  whether  the 
six  bi-monthly  correction  fields  perform  better  than  the  two 
semi-annual  correction  fields.  In  some  of  the  correction 
fields,  the  east-west  variation  is  not  large.  It  is  possi- 
ble that  a  correction  dependent  only  on  latitude  may  perform 
as  well.  Such  a  correction  would  be  easier  to  apply.  The 
numerical  ocean  prediction  experiments  necessary  to  demon- 
strate the  usefulness  of  the  corrections  are  in  progress, 
and  will    be  reported   separately. 
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Fig.  1  The  difference  (^ H-QT|  and  the  change  in  heat 
content  relative  to  200  m  (AH)  and  the  surface  heat- 
ing Q-  )  between  July  19/6  and  May  1976.  Cross- 
hatched  areas  have  insufficient  data  for  the  AH 
analysis.  Positive  values  of  ^H  indicate  increasing 
heat  content,  while  positive  QT  indicates  a  net  down- 
ward surface  heat  flux.  Units:  10*  cai  cm"z  oer  2 
months . 

-2    -J 

Fig.       2      Correction    field    (cai    cm      h   |       to    be    applied 

to  FNOC  heat  flax  fields  based  on  all  36  of*  the 
monthly  sets  of  j(AH-Qr)  evaluated  over  two-month  in- 
tervals. Negative  values  indicate  that  the  upward 
heat   flux    is   to    be    reduced   by    the   amount    shown. 


Fig,  3  As  in  Fig 
fields  for  November 
through    October    (bo 


"-g 


2,      except      seDara' 
through    April      (topi 
torn)  . 


e  correction 
and  for  May 


Fig.      4      A  s    m    i.  j.«j  .      c  , 
fields    for:    February  - 
die);    June    -   July    (bottom) 


except  bi-monthly   correcti.cn 


February  -'March    "(top)  ;    April"-    May    (mid- 


Fig.  4  (continued)  As  in  Fig.  2,  exceDt  bi-monthly 
corrections  for:  August  -  September  (top):  October  - 
November     (middle) ;     December    -    January    (bottom) . 


Fiq.    5   Time   series 

(dashed),   uncorrected  surface   neatim 

corrected  surface  heatina  (dotted!  at 


% 


heat  content 
(solid),  and 
S*N,    175*ff.    C 


is    the    corrected   total    surrace    heating.       Units    are 
cal   cm~*    over    two-mcnth    intervals. 
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Fig.  1  The  difference  (& H-QT>  and  the  change  in  heat 
content  relative  to  20O  ai  (A.3)  an^  'the  surrace  heat- 
ing (QT  )  between  July  197b  and  May  1976.  Cross- 
haxche'd  areas  have  insufficient  data  for  -he  AH 
analysis.  Positive  values  of  ^H  indicate  increasing 
hear  content,  while  oositive  QT  indicates  a  net  down- 
ward surface  heat  flux.  Units:  10%  ral  cm*1  per  2 
months . 
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Fiq.  2  Correction  field  (cal  cm  h  )  to  be  applied 
to  FNOC  heat  flux  fields  based  on  all  36  of  the 
monthly  sets  of  j(A  H-Qr )  evaluate!  ov=r  two-^onrh  in- 
tervals. Negative  values  indicate  that  the  upward 
heat  flux  is  to  be  reduced  by  the  amount  shown. 
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Pig.  3  As  in  Fig.  2,  sxrapt  separate  correction 
fields  for  November  through  April  (topi  and  for  May 
through  October  (bottom). 
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£i2«  J.  (continued)  as  m  Fig.  2,  exceDt  bi-monthly 
corrections,  tor:  August  -  Saptember  (top):  October- 
November     (middle);     December    - January  l  (bottom) . 
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Appendix   A 
MONTHLY    SETS     OF    &E-Q^,     ^H    AND    QT 

The  monthly  sets  of  &'A-QT,  &  H  and  QT  based  on  two- 
monthly  intervals  are  included,  except  for  the  May  -  July 
1976  set  which  is  in  Fig.  1.  The  form  of  these  diagrams  is 
described  in  the  caption  of  Fig.  1.  Positive  values  (solid 
lines)  are  increasing  heat  content  in  time  and  net  surface 
heating,  and  conversely  far  negative  values  (hashed  lines). 
Zero  lines  are   enhanced. 
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